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Abstract 

Density functional perturbation theory is used to analyze electron-phonon interaction in bilayer 
graphene. The results show that phonon scattering in bilayer graphene bears more resemblance 
with bulk graphite than monolayer graphene. In particular, electron-phonon scattering in the 
lowest conduction band is dominated by six lowest (acoustic and acoustic-like) phonon branches 
with only minor contributions from optical modes. The total scattering rate at low/moderate 
electron energies can be described by a simple two-phonon model in the deformation potential 
approximation with effective constants D ac ~ 15 eV and D op ~ 2.8 x 10 8 eV/cm for acoustic 
and optical phonons, respectively. With much enhanced acoustic phonon scattering, the low field 
mobility of bilayer graphene is expected to be significantly smaller than that of monolayer graphene. 

PACS numbers: 72.80.Vp, 71.38.-k, 71.15.Mb, 63.22.Rc 
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Graphene, a two-dimensional sheet of carbon atoms in a honeycomb lattice, has received 
wide attention due to its unique properties In addition to the significant interest in fun- 
damental physics, which stems in part from the relativistic-like behavior of charge carriers, 
this material is considered very promising in many applications. Particularly, the possibility 
of manipulating the band gap in a bilayer form [2] offers an additional control for nonlinear 
functionality. 

While there has been a large number of reports (both experimental and theoretical) on 
graphene and its derivatives 3|-[6), most of the studies on electron transport properties have 
concentrated on monolayer graphene (MLG). Among those requiring further investigation 
in bilayer graphene (BLG), intrinsic carrier-phonon scattering is one of the most crucial as it 
determines the ultimate limit of any electronic applications. Due to the nonpolar nature of 
the material, the deformation potential approximation is commonly used. However, the de- 
formation potential constant that quantifies the strength of electron-phonon coupling, must 
be determined outside the developed theory, in the form of an empirical fitting parameter. 
The experimental results lie in a broad range (e.g., D = 10 — 50 eV) 4], 6] and further 
improvement of accuracy is essential for reliable analysis of transport characteristics. The 
ab initio numerical methods provide an alternative approach in addressing this problem as 
demonstrated successfully in MLG (3, O]- In this study, we conduct a first principles cal- 
culation of electron-phonon interaction in BLG with AB (also known as Bernal) stacking. 
Specifically, the density functional perturbation theory (DFPT) is used to obtain the crystal 
phonon spectra and electron-phonon coupling matrix elements. The results illustrate the 
differences with MLG, indicating significant discrepancies in transport properties. 

Although it has been reported that the electron-phonon matrix elements in MLG can be 
estimated with an excellent accuracy using the GW approximation [?| for the phonon modes 
associated with the observed Kohn anomalies jsj], a considerably longer (compared to the 
DFPT) calculation time that the GW method requires, as well as the fact that the accuracy 
of DFPT is sufficient in many applications including carrier transport, justifies the use of 
DFPT for calculating electron-phonon matrix elements in BLG. We also demonstrate that 
the Kohn anomalies that lead to the largest discrepancy between DFPT and GW results, do 
not play an important role in the intrinsic electron scattering in BLG and so the accuracy 
of the calculated matrix elements is not essential in this case. 

Figure [T] shows the calculated phonon dispersion along the T-K direction. Since 
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there are four carbon atoms in the unit cell, BLG has 12 phonon branches. However, only 
six lines are distinguishable across the entire first Brillouin zone (FBZ) as each of them 
are doubly degenerate due to the weak van der Waals coupling between the two carbon 
layers. The only exception, when the splitting has an appreciable magnitude, is the case 
of ZA (acoustic out-of-plane) and ZO' ("layer breathing") modes near the FBZ center. In 
addition, two other modes deviate from the pure acoustic nature with very small non-zero 
frequencies at the T point as shown in Fig. QJb) (denoted as TA2 and LA2). Hence, only 
three branches are truly acoustic (ZA, TAl, LAI), while we will use a term "antisymmetric 
acoustic" to refer to the three low-frequency counterparts with acoustic-like behavior (ZO', 
TA2, LA2). Their polarization eigenvectors are akin to those of the acoustic modes in MLG, 
with the atomic displacements in two layers having opposite phase (hence "antisymmetric"). 
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A similar terminology has been used in Raman spectroscopy studies 

The electron spectrum of BLG has two close conduction bands (7r* and tt^); the 
distance between the bottoms of 7r* and 7r| is found to be in a range 0.35 — 0.4 eV 



121 ] . The accurate estimate of the intrinsic electron scattering in BLG requires taking into 



account all four possible scattering scenarios: intraband (nl — > irl, tt^ — > ^2)1 an d interband 
(-7T* T^- ^2). Correspondingly, we calculate four sets of matrix elements | fi'k+q^k I ^ na ^ represent 
the quantum mechanical probability of electron transition from the state (i, k) to (j, k + q), 
where i and j are the initial and final energy band respectively. The matrix element is 
calculated only for wave vectors q in the irreducible Brillouin zone. However it is easy to 
obtain the matrix elements in the entire FBZ, by restoring the star of each q- vector from the 
irreducible zone (which requires applying corresponding symmetry operations to the vectors 
q and k simultaneously). The result of this calculation for intraband scattering (when the 
electron is at the bottom of 7r^ , i.e. k = K) is shown in Fig. |2j It is not surprising that MLG- 
like phonon modes (ZA, TAl, LAI, ZOl, TOl, LOl) and their BLG-counterparts (ZO', TA2, 
LA2, Z02, T02, L02) have very similar electron-phonon matrix elements throughout FBZ. 
The difference is only exhibited in some cases, in points of high symmetry (q = K and 
q = 0). This can be explained by selection rules that apply to phonon- mediated electron 
transitions, and it results from different symmetries of atomic displacements of these two 
sets of modes ljij]. As compared to MLG, the matrix elements have the same order of 
magnitude. The essential difference between matrix elements of MLG and the ones shown 
in Fig. [2] is that none of the in-plane optical phonon modes of BLG reveal Kohn anomalies 
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that were reported to be observed in a form of singularities at points q = K and q = jg]. 
However, further calculations have shown that strong peaks similar to those found in MLG 
occur in the corresponding matrix elements of interband scattering 7r* — > 7r|. Figure E] shows 
the matrix element of interband scattering 7r* — > H2 due to the phonons of branches 9 and 
10 (TOl and T02). The symmetry analysis of these two branches imposes a selection rule 
that allows only participation of one of the two phonons with q = K [13] . In agreement 
with this prediction, the sharp peaks at q = K are observed only in the case of TOl. At 
the same time, the other Kohn anomaly (at q = 0) is present in matrix elements of all four 
MLG-like optical branches: TOl, T02, LOl, L02 (see Fig.E}. 

Based on the obtained matrix elements, the electron scattering rates are calculated using 
Fermi's golden rule 

;) k =|E|«£&| <f. W 

where A { ^ u = (JV£ + § ± |) 6 (£j +q - E l k ± huj£) , N£ is the phonon population factor, 
and the sign plus (minus) corresponds to the emission (absorption) of the phonon oj^. The 
temperature dependence of the scattering rate stems only from the the numbers of phonons 
as we sum over wave vectors q in the FBZ. Preliminary estimates confirm that the scattering 
rate, as a function of the initial electron state k in the vicinity of the Dirac point K, exhibits 
isotropy. Therefore we only consider the K — V direction in k-space and can represent r _1 
as a function of electron energy. 

In the integration over FBZ the tight binding approximation was used to describe the 
electron energy spectrum in BLG The second conduction band (ir^) has nearly the 
same shape as the lowest conduction band (7r*) near Dirac point K: it is parabolic in close 
vicinity of k = K and quasi-linear beyond that. The quasi-linear approximation holds true 
for quite high electron energies and therefore, this description appears sufficient in most 
calculations. In the tight binding approximation, the distance between the bottoms of 7rJ 
and 7i~2 is 71, which is the interlayer hopping integral in BLG (we use standard notations 



known as the Slonczewski-Weiss-McClure parametrization model 14|,|l5|]). A recent resonant 



Raman studies show the result 71 ~ 0.4 eV 
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12j , which is the value used in our calculation. 



Figure H] shows the total scattering rates associated with all four possible scenarios in- 
volving two conduction bands. Since the recent Raman experiments in BLG have shown 
the Kohn anomalies located in the q-points of high symmetry, one might expect their sig- 
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nificance in the overall picture of electron-phonon scattering in BLG. The separate analysis 
of different phonon branches in BLG, however, shows a qualitatively different picture. In 
case of intraband scattering in the lowest conduction band (71"*) the scattering rate due to 
absorption or emission of MLG-like optical phonons (ZOl, Z02, TOl, T02, LOl, L02) is 
negligibly small in a wide range of energies (E < 0.7 eV). An explanation for this result 
comes directly from comparison of corresponding matrix elements of MLG [8j] and those of 
BLG for irl — > -k\ scattering: in the latter case the sharp peaks are missing. The Kohn 
anomalies similar to those in matrix elements of MLG, were observed in BLG - in the in- 
terband scattering — > 7r| (see Fig. [3]). However, the corresponding scattering rates have 
proven to be negligibly low. This result can be understood if one compares electron dis- 
persions in MLG (71"*) and BLG (7r* and ir^). A simple analysis demonstrates that electron 
transitions 7r* — > ~k\ involving phonons q = and q = K (where the anomalies occur) are 
prohibited by the energy and momentum conservation laws. In other words, the contour 
of integration from Eq. (TjQ) never crosses (or comes close to) the peaks at the points T and 
K observed in the matrix element of BLG. The estimated minimal distance between this 
curve (which can be approximated by a circle) and the point q = K, is 0.034 (2ir/a) for 
emission and 0.014 (2n/a) for absorption. For comparison, the analogous distance in case 
of intervalley scattering due to TO phonons in MLG, is much smaller, ~ 10 -7 (27r/a), which 
ensures the dominant role of the corresponding modes in electron scattering. The same con- 
sideration shows that the TOl phonon with q = cannot influence the interband scattering 
in BLG. As a result, the total scattering rate in BLG, when the electron is initially in ttI 
(i.e., 7r^— n^}) is somewhat smaller at higher energies (E^ > 200 meV), if compared to 
the total rate in MLG (see Fig. H]). 

On the other hand, the low-energy scattering rate in BLG is substantially higher, which 
can be explained by a larger density of states in the vicinity of the Dirac point. More 
specifically, the density of state tends to zero linearly for the massless Dirac fermions of 
MLG while it is a constant in BLG as k — > K. This non-zero density of states at the 
bottom of ttI and tt^ leads to a peculiar effect. In MLG, the intravalley scattering on 



acoustic phonons can be treated as quasielastic as it was shown previously |16l . |17| ; only 
the modes with q — > can participate. As a result, the scattering rate is a linear function 
of energy that vanishes as 0. In BLG, however, the intravalley scattering with 

absorption of acoustic phonon shows a different behavior. Aside from q = 0, the energy and 
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momentum conservation in the electron scattering allows an additional solution - a circle 
with radius |q| = q c = Vg^i/hvp. In the case of antisymmetric acoustic modes, an equivalent 
solution exists - a circle of radius q c = y jxour / fiv^, where ur is the non-zero frequency of 
an antisymmetric acoustic mode in T-point. This leads to a non-zero absorption rate even 
at zero electron energies. It is important to emphasize that even though the deformation 
potential approximation can be used to describe electron-phonon coupling in BLG [6), the 
intravalley scattering cannot be treated as quasi-elastic due to the additional solution with 
q^O. The estimated value of the deformation potential constant in BLG is D 15 eV, if 



we chose the TA2 branch as dominant 



181 ]. This is much closer to the value obtained for 



graphite (D « 16 eV Jl9|) than that for MLG (D ^4.5-5 eV [8]). The scattering rates 
in BLG and graphite j2o| turn out to be quite similar as well. The same dominant intrinsic 
scattering mechanism in graphite and BLG is intravalley transitions due to absorption and 
emission of acoustic phonons, while the contribution of the optical phonons is of secondary 
importance. In particular, D op rj 2.8 x 10 8 eV/cm for the optical phonons in BLG, which is 



almost one order of magnitude smaller than in MLG [8J, Il8j. With much enhanced acoustic 
phonon scattering at low energies, the low field mobility of bilayer graphene is expected to 
be significantly smaller than that of monolayer graphene. 

In summary, we have studied the electron-phonon interaction in BLG with AB stacking 
and found that the scattering rates exhibit qualitatively different behavior, if compared 
to MLG. As a result of a different electron spectrum (and thus, the density of states) in 
conduction bands of BLG, the intraband scattering rates due to absorption of acoustic 
and antisymmetric acoustic phonons at low electron energies have quite large values (~ 
10 9 — 10 10 s _1 ). The calculations confirmed that this scattering mechanism remains dominant 
in BLG, at all relevant electron energies in ir^ — the property that is not observed in MLG. 
In 7i"2, the scattering rate is approximately a constant, ~ 5 x 10 12 s _1 . The overall picture of 
electron-phonon coupling in BLG has more in common with that of the bulk graphite than 
it does with MLG. The estimated values of the deformation potential constants are D ~ 15 
eV and D op 2.8 x 10 8 eV/cm, for acoustic and optical phonons respectively. 
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FIG. 1: (Color online) The phonon spectrum of BLG along the line T — K, obtained by DFPT 
method (a). There are total of 12 phonon branches in BLG - twice as many as in MLG. But 
the weak van der Waals interlayer coupling results in effective double degeneracy of every branch, 
everywhere in FBZ except the vicinity of T-point where a small splitting is exhibited in case of 
first six branches (b). The largest splitting occurs between two out-of-plane branches: ZA and ZO' 
("breathing mode"). 
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FIG. 2: (Color online) Intraband (ir*) matrix elements 
DFPT for k at the conduction band minimum (i.e., the Dirac point) as a function of phonon 
wavevector q and branch number v. If compared to the matrix elements of MLG, the main 
difference is the absence of Kohn anomalies - sharp peaks in q = K (TOl, T02) and q = (LOl, 
L02). 
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FIG. 3: Interband matrix elements (tt^ — > 7r|) (in units of eV) for phonon branches TOl and T02 
(k = K). The Kohn anomalies are clearly visible in a form of sharp peaks. The observed picture 
agrees with the selection rules for electron-phonon interaction in BLG. [2] 
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FIG. 4: (Color online) Scattering rates at T = 300 K as a function of electron energy E^, in BLG. 
The results are rather different depending on whether the electron is initially in the conduction 
band 7r* (a), or 7r| (b). In the former case, the dominance of the intraband (7r*) scattering is 
clear. The scenario when the electron is initially in ir^ becomes more relevant when the gate bias 
is applied. 
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